Full Current Statistics in the Regime of Weak Coulomb Interaction 
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We evaluate the full statistics of the current via a Coulomb island that is strongly coupled to the 
leads. This strong coupling weakens Coulomb interaction. We show that in this case the effects 
of the interaction can be incorporated into the renormalization of transmission eigenvalues of the 
scatterers that connect the island and the leads. We evaluate the Coulomb blockade gap in the 
current-voltage characteristics, the value of the gap being exponentially suppressed as compared to 
the classical charging energy of the island. 
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The electric charge is quantized. There are two impor- 
tant manifestations of this fact in electron transport on 
mesoscopic scale. 

The first manifestation is a shot noise (ij. The quan- 
tum mechanics modifies the shot noise in an arbitrary 
mesoscopic conductor with respect to its classical value 
S = 2el. It also changes the full statistics of charge 
transfer via a mesoscopic scatter, so that the statistics is 
not the simple Poissonian one 

The second manifestation of electron charge quan- 
tization is the Coulomb blockade phenomenon. The 
Coulomb blockade is most strong provided the conduc- 
tance G of a mesoscopic system is much smaller than the 
conductance quantum Gq = e^/27r?i 0. The statistics 
of charge transfer in this strongly interacting case is that 
of a classical stochastic Markov process 0] . On the con- 
trary, quantum mechanics is of importance for Coulomb 
blockade in the opposite limit G ^ Gq where Coulomb 
interaction is weakened. 

The charge quantization manifests itself in both 
ways in full current statistics under conditions of weak 
Coulomb interaction, this provides the motivation for the 
present study. 

It has been understood that the charge quantization 
survives even in the limit G > Gq H H Hll US • The 
ground state energy of the Coulomb island was shown to 
retain the periodic dependence on the induced "off-set" 
charge q, thus indicating the Coulomb blockade. How- 
ever, due to large quantum fluctuations of charge in the 
island, the effective charging energy Ec (defined as the 
periodic q-dependent part of free energy) is exponentially 
suppressed as compared to the classical charging energy 
Ec = e^/2G. It is important that the weak charge quan- 
tization persists not only for tunnel junctions 0, 0] 
connecting the island and the leads. The connectors can 
be arbitrary mesoscopic scatterrers The quantiza- 

tion completely vanishes only for constrictions with per- 
fectly transmitting channels H, 0| . 

This research was concentrated on the equilibrium 
properties, an alternative is to study transport. Recent 
studies link the short noise in the conductor to the inter- 



action correction to the conductance. The interac- 
tion correction to full current statistics was analyzed in 
[l^ for a scatterer embedded in the electromagnetic envi- 
ronment. The relation between interaction correction to 
the conductance and the formation of Coulomb blockade 
in an island was investigated in |l3l |. 

In this Letter we evaluate the full statistics of the cur- 
rent via the Coulomb island defined by several meso- 
scopic scatterers assuming weak Coulomb interaction, 
G ^ Gq. The results can be summarized as follows. 
At energy scale E ^ goEc, go = Gq/Gq being the di- 
mensionless conductance of the system in the absence of 
interaction, the dominant effect of Coulomb interaction 
is the energy-dependent renormalization of the transmis- 
sion eigenvalues TA of the mesoscopic scatterers labeled 
by fc. 
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The renormalization of a similar form was previously ob- 
tained in [l8| for a scatteter in the weakly interacting 
Id gas and in for a single multi-channel scatterer 
shunted by an external impedance. We thus prove this 
simple relation for a Coulomb island. The full current 
statistics is readily obtained from the energy-dependent 
Tn with using non-interacting scattering theory ap- 
proach of 0, E3- This gives the voltage dependence 
of conductance, shot noise and all higher cumulants of 
charge transfer. In contrast to the case of a single scat- 
terer, the renormalization of all transmission eigenval- 
ues may break down at finite energy Eq oc goEce~"^° , 
a being a numerical factor depending on the details of 
the initial transmission distribution. Remarkably, Ec 
coincides with the equilibrium effective charging energy 
evaluated with instanton technique. 10] However, the 
renormalization stops at the effective Thouless energy 
E-Th ~ G{E)S/Gq, 5 being mean level spacing in the 
island. This gives rise to two distinct scenarios at low 
enegy. If go > a^^ \n{Ec/S), Coulomb blockade does 
not occur with zero-bias conductance being saturated at 



2 




\—0X 



FIG. 1; Multi-terminal (M=3) Coulomb island. 



the value G(£'Th) > Gq. Alternatively, G{0) w and 
Eq defines the Coulomb gap. 

Let us give the details of the model in use. The 
Coulomb island is characterized by two parameters: 
charging energy Ec and the mean level spacing S, Eq 3> 
6. The island is connected to Af > 2 external leads by 
means of M arbitrary mesoscopic scatterers. (Fig. 1) 
Each scatterer, or connector, i is characterized by the set 
of transmission eigenvalues Tn^. We assume that the is- 
land is strongly coupled to the leads, go = J2n m ^ ^■ 
Our goal is to evaluate the cumulant generating func- 
tion (CGF) S{[xi]))- The Fourier transform of exp(— S") 
with respect to counting fields Xi gives the probability 
P{{Ni}) for Ni 1 electrons to be transferred to the 
terminal i during time interval to. (See 0) The deriva- 
tives of S give the average value of currents, shot-noise 
correlations and higher order moments of charge transfer. 

To evaluate the CGF of the Coulomb island, we have 
extended the semiclassical approach for the FCS of the 
non-interacting electrons [l5j- To account for Coulomb 
interactions, we introduce a dynamical phase variable 
(p{t) that results from the Hubbard-Stratonovich 
transform of the charging energy. Its time derivative, 
(p{t)/e, presents the fluctuating electrostatic potential of 
the island. The CGF S{{xi}) can be then represented 
in the form of a real-time path integral over the fields 
fi,2{t) residing at two branches of the Keldysh contour 



+ OC 



exp 



{~S{X^}) = J D^,,2{t)cxp['-Ec' J dt{^\-^l) 

— oo 

^5[i({G,G^}) -*^5-iTr{(*a, - <i>)G}} (2) 



Here $ 



is the matrix in Keldysh 
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space, 2x2 matrix G(ii,i2) presents the electron 
Green function in the island that implicitly depends on 
<yJi,2(i)- The trace operation includes the summation over 
Keldysh indices and the integration in time. The contri- 
bution of each connector S'ioL has a form found in the 



circuit theory [isl 
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{G, G^} denoting the anticommutator of the Green func- 
tions with respect to both Keldysh and time indices. 
The Green functions in the leads Gfc(Y) are obtained by 
X-dependent gauge transformation [ig of the equilib- 
rium Green functions in the reservoir fc, Ct-^\ G^(e) = 



exp(«XfcT3/2)G[.°^(e) exp(-jxfcr3/2), where g{°' are given 
1-2/fe -2/fc 
-2(1 -/fc) 2fk-l 
the electron distribution function in the fc-th reservoir. 
The expression ^ is valid under assumption of instanta- 
neous electron transfer via a connector, this corresponds 
to energy independent t]^^ . 

In order to find G(ii,i2) at given (pi, 2(^)1 we minimize 
the action with respect to all G{t 1,12) subject to the con- 
strain G o G — 5{ti — 12). This yields the saddle point 
equation for G(ti,t2): 

f''[G^,G] 
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where [..,..] denotes the commutator in the Keldysh- 
time space. This relation expresses G(ii,t2) = 
G{ti,t2; [<j5i_2(0]) via the reservoir Green functions G^'^l. 
This circuit theory relation is similar to obtained in |l5l| . 
It disregards the mesoscopic fluctuations, since those lead 
to corrections of the order of ~ 1/go at all energies, 
whereas the interaction corrections are of the order of 
~ l/(?oln(i?) tending to diverge at small energies. If 
fi.2{t) = 0, Eq. Q separates in energy representation 
and coincides with that of Ref. |3 . 

This sets the model. We start the analysis of the 
model with perturbation theory in (/?i_2 around the semi- 
classical saddle point G{ti,t2) ~ Go, >fi,2{t) = 0. The 
phase fluctuations are small, <5(/3^ ~ 1/go, so we keep only 
quadratic terms to the action The resulting Gaus- 
sian path integral over ipi^2 can be readily done. This 
procedure is equivalent to the summation of all one-loop 
diagrams of the conventional perturbation theory, i.e. to 
the "random-phase approximation" (RPA). 

For the rest, we restrict ourselves to the most inter- 
esting low voltage/temperature limit, maxje^, fcT} <C 
goEc- In this limit, we evaluate the interaction correc- 
tion to the CGF with the logarithmic accuracy. It reads 

A^, = ^lnf 'f,^,^, )^ (5) 
^ 50 \max{eV, kT} J ^ ' 
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n, k 
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provided max{ey, fcT} > -Exh, where i?Th = ffo^ is the 
Thouless energy of the island. In the opposite case. 
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max{eV, kT} < £^Th, the voltage/temperature should be 
replaced with i?Th- Note, that the correction ((SJ is con- 
tributed by only virtual inelastic processes that change 
the probabilities of real elastic scatterings. 

For simplicity, we consider the shot-noise limit eV ^ 
kT only. Then the magnitude of the correction shall be 
compared with the zero-order CGF ~ toeVgo- This 
implies that the perturbative RPA result ((S} is applicable 
only if gQ^\n{goEc/eV) <C 1. At lower voltages AS^ 
logarithmically diverges. This indicates that we should 
proceed with a renormalization group (RG) analysis. 

We perform the RG analysis of the action (O in the 
one-loop approximation, this corresponds to the "poor's 
man scaling". This is justified hy g ^ 1. We follow 
the conventional procedure and decompose ip{t) onto 
the fast (fif and slow parts (ps- On each step of RG pro- 
cedure we eliminate the fast degrees of freedom in the 
energy range E — 5E < uj < E to obtain new action 
Se-se[^s], E being the current ultraviolet cutoff. Our 
key result is that the change in the action at each step of 
RG procedure can be presented as a change of transmis- 
sion eigenvalues Tn^ . 

Therefore, the RG equations can be written directly 
for transmission eigenvalues and take a simple form 
The equations are to be solved with initial conditions 
at the upper cutoff energy E = goEc, those are given 
by "bare" transmission eigenvalues Tlt\E = g^Ec) = 

\k] 

Tn ■ The RG equations resemble those for the trans- 
mission coefficient for a scatterer in the weakly inter- 
acting one-dimensional electron gas |0| and for a single 
multi-channel scatterer in the electromagnetic environ- 
ment Il2j. The effective impedance Z is just replaced 
by inverse conductance of the island to all reservoirs, 
G{E) = GQY,n,kTn\E). The important difference is 
that this conductance is itself subject to renormalization. 
The difference becomes most evident in the case when all 
contacts are tunnel junctions, <C 1. In this case, one 
can sum up over k,n in Eqs. JQ) to obtain the RG for 
the conductance only : dG/dhiE = 2Gq. This renor- 
malization law [ioj was recently applied to conductance 
of granular metals. The Eqs. could be also derived in 
the framework of functional RG approach to cr-model of 
disordered metal. [iH . 

We solve the RG Eqs. ^ in general case to obtain 

tV^He) = ll^-ljy/(l-rW(l-y)), (6) 

Hg,Ec/E) = -\Y.Hl~T!^\l-y)) (7) 

n, k 

The first equation gives the renormalized transmission 
eigenvalues at a given value E of the upper cutoff in 
terms of variable y{E), < ?/ < 1. The second equa- 
tion implicitly expresses y{E). 

We note that the energy dependence of transmission 
coefficients induced by interaction is very weak provided 
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FIG. 2; The total conductance of the Coulomb island versus 
the energy: two scenarios. We assume ln(_Bc/<5) = 10.0. Ar- 
rows show the energy scale ~ 5. Pane (a): tunnel connectors, 
go changes from 42 (upper curve) to 14 (lowermost curve) 
with the step 4. Pane (b): diffusive connectors, go changes 
from 18 to 6 with the step 2. 

G{E) ^ Gq: If energy is changed by a factor of two, the 
conductance is changed by ~ Gq. To use the equations 
for evaluation of FCS at given voltages V^'^^ of the leads, 
one takes T]t\E) at upper cutoff E = max^ (yW), and 
further disregards their energy dependence. Then one 
can follow the lines of Ref. [l5j: It is convenient to in- 
troduce the function 5'W(a;) = - E„ + \Tn\x - 1)] 
to incorporate all required information about transmis- 
sion eigenvalues. The renormalization of S^^^ in terms 
of y is especially simple: S'^^^x.y) = S^''\{x + l)y — 
1) - 5^(2?/ - 1). From this one readily finds the con- 
ductance of each scatterer, G^'^^y) = 2GQdS^''^ /dx{l,y), 
as well as the renormalized transmission distribution 
T^p^^\T,y) = (2/7r)Ini{9S'W/9a;(l-2/r-iO,?;)}. 

The RG equations (1) have a fixed point at tIi'^ ~ 
0, y = that occur at finite energy 

E^Ec^g,Ec\{{l-TV^^f'^ (8) 

k,n 

This indicates the breakdown of RG and formation^ of 
Coulomb blockade with the exponentially small gap Ec- 
The same energy scale was obtained from equilibrium 
instanton calculation of Ref. 10]. For a field theory, one 
generally expects different physics and different energy 
scales for instantons and perturbative RG. The fact that 
these scales are the same shows a hidden symmetry of 
the model which is yet to understand. 

Alternative low-energy behavior is realized if the cur- 
rent cut-ofl[ reaches E'xh = G{E)S/Gq. (Fig.2) The log 
renormalization of the transmission eigenvalues stops at 
this point and their values saturate. We thus predict 
a sharp crossover between the two alternative scenar- 
ios, that occur at value of go — gc corresponding to 
Ec — This value equals gc = a~^\n{Ec/5), where 

~ X^n k ln(l — tI'''), and depends on transmission 
distribution of all connectors. If all connectors are tunnel 
juctions, ut = 2. For diffusive connectors, an = 7r^/8 
and the energy dependence of the total conductivity is 
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land. We have revealed that the interaction effect can be 
incorporated into a simple energy-dependent renormal- 
ization of transmission eigenvalues, this enables easy eval- 
uation of all transport properties. The Coulomb block- 
ade develops only if the total conductance of the island 
is below some critical value ^ GQ\og{Ec/S), otherwise 
the interaction correction to transport saturates at low 
energies. 

This work is a part of the research program of 
the "Stichting voor Fundamenteel Onderzoek der Ma- 
terie" (FOM) , and we acknowledge the financial support 
from the "Nederlandse Organisatie voor Wetenschap- 
pelijk Onderzoek" (NWO). 



FIG. 3: An example of FCS: Probability to measure equal 
currents in the 3-terminal island with identical junctions. The 
curves (a-d) present diffusive connectors at different values of 
the renormalization parameter y{eV): (a) - 1.0, (b) - 0.5, (c) 
- 0.1, (d) - 0. No interaction effect is seen for tunnel (1) and 
ballistic (2) connectors. 



given by goiV) - ffov^ctgv^, C = 2gg^ h-i{goEc /eV). 
(Fig. 2) 

To give an example of FCS calculation in the regime 
of weak Coulomb interaction, we consider the 3-terminal 
Coulomb island with identical junctions, discussed in [isj 
(Fig. 1). The voltages apphed are V3 = V, Vi,2 = We 
plot in Fig. 3 the log(P), the logarithm of the probability 
to measure the same currents to the terminals 1 and 2, 
Ii — I2 — —1^/2, versus the current I3 measured in the 
terminal 3 (Fig. 3). Both log(P) and I3 are normalized 
by the average current < J3 >, so that in the absence of 
interaction the curves corresponding to different voltages 
are the same (assuming eV T), the shape of the curve 
reflecting the transmission distribution in the contacts. 
To take interaction into account, we change the trans- 
mission eigenvalues according to Eqs. (6,7) and evaluate 
the probability with the method of Ref.|l5|. 

The curves 1 (tunnel junctions) and 2 (ballistic con- 
tacts) stay the same not depending on the renormaliza- 
tion. Indeed, the renormalization does not affect ballistic 
transmission, as to tunnel junctions, it affects only their 
conductances. The interaction effect is visible for diffu- 
sive junctions. The curves (a)-(d) correspond to different 
values of y{E = eV). The transmission distribution of 
each contact evolves from the diffusive form {po{T) = 
g{y)/2TVl - T) at the highest voltage (y « 1) to the 
double junction form (pdjCT) = g{y)/-KT^^'^Vl^T) at 
the lowest voltage (y « 0) This visibly changes the 
form of logP(/3). 

To conclude, we have analyzed the effect of weak 
Coulomb interaction (G ^ Gq) on FCS in Coulomb is- 
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